function CrackZone
global glb
Sy=glb.mater.Sy;
KI=glb.crack.KI;
glb.crack.R=1/2/pi*(KI/Sy).^2;